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Abstract 

We present the whole set of equations with regularity and match- 
ing conditions required for the description of physically meaningful 
stationary cylindrically symmmetric distributions of matter, smoothly 
matched to Lewis vacuum spacetime. A specific example is given. The 
electric and magnetic parts of the Weyl tensor are calculated, and it 
is shown that purely electric solutions are necessarily static. Then, 
it is shown that no conformally flat stationary cylindrical fluid exits, 
satisfying regularity and matching conditions. 
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1 Introduction 



In a recent paper [1] we have presented a systematic study of static cylin- 
drically symmetric matter distributions. To justify such endeavoiu we men- 
tioned, on the one hand the still ongoing discussion about the precise meaning 
of the independent parameters of the corresponding vacuum spacetime, and 
on the other, the renewed interest in cylindrically symmetric sources in rela- 
tion with different, classical and quantum, aspects of gravitation (see [2] and 
references therein). 

It is clear that the extension of the previous study to the stationary case, is 
not only justified by the same arguments, but in addition allows for studying 
the role of the rotation of the source. 

However, it should be stressed that in spite of the simplifications in- 
troduced by the cylindrical symmetry, a purely analytical handling of the 
problem requires further restrictions. Accordingly, we shall consider here the 
shear-free (rigid rotation) case. 

As in [1] we want to bring out the relationship between the Weyl tensor 
and different aspects of the source. This last question is in turn motivated by 
the very conspicuous link existing in the spherically symmetric case between 
the Weyl tensor, the inhomogeneity of the energy density and the anisotropy 
of pressure [3]. For doing so we have calculated the Weyl components as 
well as the components of its electric and magnetic parts. Unfortunately, no 
such link can be established in our case. However it is shown that purely 
electric solutions are necessarily static. This last result, together with results 
obtained in [1] imply that no conformally flat stationary cylindrical fluid 
exists, satisfying regularity and matching condition. 

The paper is organized as follows: in section 2 we present the general form 
of the energy momentum tensor, the line element, the Einstein equations, the 
regularity conditions and the active gravitational mass. The exterior space- 
time as well as junction conditions are discussed in section 3. In section 4 a 
speciflc example, which we believe is new, is presented. In section 5 the elec- 
tric and magnetic parts of Weyl tensor are presented and the non existence 
of conformally flat models satisfying Darmois conditions is established. Fi- 
nally, some conclusions are presented in the last section and the components 
of Weyl tensor are given in an appendix. 
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2 Inside the source 



We shall first describe the kind of source we are deahng with as well as its 

spacetime. Then we shall write the corresponding Einstein equations and dis- 
cuss the regularity conditions to assure the correct behaviour of the solutions. 
Finally we present a very compact expression for the active gravitational mass 
of the source. 

2.1 The source 

We consider a stationary cylindrically symmetric anisotropic non-dissipative 
fiuid bounded by a cylindrical surface E and with energy momentum tensor 
given by 



Tc^P = (/X + PrWaVp + PrQaP + (^0 ' Pr)K^Kp + (P, - Pr)S^Sp, (1) 



where, ji is the energy density, P^, Pz and P^ are the principal stresses and 
Vai Ka and Sa are four- vectors satisfying 



2.2 The interior spacetime 

We assume for interior metric to E the general stationary cylindrically sym- 
metric line element, which can be written 



where /, A;, 7 and / are all functions of r. To represent cylindrical symmetry, 
we impose the following ranges on the coordinates 



We number the coordinates = t, = r, x"^ = z and x^ = (f) and we choose 
the fluid being at rest in this coordinate system, hence from (2) and (3) we 
have 



(is2 = _fdf + 2kdtd(j) + e^{dr^ + dz'^) + ld<f, 



(3) 



—00 < i < 00, < r, —00 < z < 00, < < 27r. 



(4) 
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It is worth noticing that we are assuming that our coordinates are co- 
rotating with the fluid {V" = 5o/v7)) which imphes, as it can be very easily 
shown, that the fluid is shear-free. Thus we shall consider a shear-free, 
expansionless fluid distribution. 



2.3 The Einstein equations 

For the Einstein field equations, G^/? — i^T^p, with (1), (3) and (5) we have 
the non-null components 

Coo = [2fDS' + ^fDD" - 2D'f + 2DD'f 

-Sfifl' + k")] = Kiif, (6) 
^ '2kD^Y + 4kDD" - 2D^k" + 2DD'k' 

-3k{f'l' + k'^)] = -Kuk, (7) 



4L>2e7 



Gu = -^.{'^DD'i + fl' + k'^) = KPre\ (8) 
G22 = -^{2DD'i - ADD" + /'/' + k'^) = KP,e\ (9) 



G: 



1 



33 



AD'^e^ 



2lD^i' + AIDD" - 2DH" + 2DD'l' 



K 



-3l{fl' + k")\ = jifik' + P^D'). (10) 

where the primes stand for differentiation with respect to r and 

D^^fl + k\ (11) 

Thus we have five equations for the eight unknown functions, namely /, 
k, 7, I, ji, Pr, Pz and P^. Therefore three equations of state relating the 

matter variables or ad hoc assumptions on the metric functions have to be 
imposed in order to integrate (6)-(10) (e.g. see the example in section 4 
below) . 

';/3 



The Bianchi identity T^.q = from (1) reduces to. 



P; + ^(/. + Pr)j + (Pr -P4>)^ + liPr " P.W = 0. (12) 



From (6) and (7) we have 

fcGoo + /Go3 = ^(^^^^)' = 0, (13) 

which allows to be integrated, producing 

kf - fk' = CiD, (14) 
where Ci is an integration constant. 

2.4 Regularity conditions 

Regularity conditions on the symmetry axis imply [4] 

hm S^^^L^ ^ 1 (15) 
r-o AX ^ ' 

with X — g^^. Then (15) yields 

The requirement that gf^^ = on the axis, gives 

I ^ 0, (17) 

where = means that the equation is evaluated at r = 0. 

Since along the axis physically there cannot be singularities, we impose 
that spacetime tends near the axis to a flat spacetime, hence we scale the 
coordinates so that as r — > we have 

ds'^ = -dt^ + 2ujr^dtd(l) + dr^ + dz^ + r'^dcf)'^, (18) 

from which it follows that 

/=^1, A;^7°0, (19) 



implying 



D = 0, (20) 
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and from (16) and (18) 

I' = 0. (21) 

Finally, from the above and the requirement that Einstein tensor components 
in (6-10) do not diverge, it follows that 



f^k'^ k" - k'^ ^ 0. (22) 
2.5 The Whittaker mass 

The Whittaker formula [5] for the active mass per unit length m of a station- 
ary distribution of perfect fluid with energy density /j, and principal stresses 
Pr, Pz and inside a cylinder of surface S is 

m = 27r / (/X + + + P^)^/^dr, (23) 

J 

where g is the determinant of the metric. Now substituting (3) and (6-10) 
we obtain 



m 



where we have used (14). Integrating (24) and using the regularity conditions 
on the axis we obtain 

^.2.flr-c^ (25) 

3 Exterior spacetime and junction conditions 

For the exterior spacetime of the cylindrical surface S, since the system is 
stationary, we take the Lewis metric [6] and consider its Weyl class [7], 

ds^ = -Fdt^ + 2Kdtd(i) + e^{dR^ + dz^) + Ld(/)^ (26) 

where 

F = aR'-'' - aS^R'+'', (27) 

X = -(1 - ab5)SR^+'' - abR^"", (28) 

= p(-'-i)/2, (29) 

L = iln^i^i+n _ ab'R'-^, (30) 



with 

5 = -, (31) 
an 

and a, b, c and n are real constants. The coordinates t, z and can be taken 
as in (3) and with the same ranges (4). The radial coordinates in (3) and 
(26), r and R, are not necessarily continuous on E in order to preserve the 
range of (f) in (4). Also, since we are considering an expansionless cylinder and 
the fluid is at rest in our coordinate system, the equation of the boundary 
surface should read r = constant and R — constant, from inside and outside, 
respectively. 

In accordance with Darmois junction conditions [8] we obtain that the 
metric coefficients (3) and (26) and its derivatives must be continuous across 

f^F, k = K, 7 = r, l = L, (32) 
/'si S^R^ + 

/ " i^ + ^'^^^i+^-i?!-' ^^^^ 



A;' E 1 (1 - abS)5R'' - abR-"" 

J'R^ '^{l-ab5)5R^+^ + abR^-^' ^ ^ 



7 ~ ^ ^(1 - ab5yR^+^ - a^m^-^' ^ ' 

hence we have from (8) Pr — 0, as expected, and for (14) and (25) 

ci = 2c, (37) 
m^a+^, (38) 

where a is the Newtonian mass per unit length, as measured outside the 
source, given by 

(39) 

The two parameters b and c are responsible for the stationarity of the space- 
time, accordingly the second term in (38) may be interpreted as the "rota- 
tional" contribution to the active gravitational mass of the source. 
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For a cylindrically symmetric rigidly rotating dust, van Stockum solution 
[9], one obtains 

*c=-ii^. (40) 

or with (38-40) 

m = a{l-2a), (41) 

which could be interpreted as if, in this case, rotation diminishes the active 
gravitational mass. At first sight this result seems counterintuitive, as we 
should expect the energy of the rotational motion to increase the total energy 
of the source. 

However, the situation is not so simple. First of all observe that the van 
Stockum solution has no static limit (see next section). Indeed, reducing 
its rotation to zero amounts to delete the source (the spacetime becomes 
flat). Therefore it is not clear that in this case, m as given by (38) might be 
interpreted as describing the static plus the "rotational" contributions to the 
active gravitational mass of the source, since no static contribution exists for 
this case. 

In the same line of arguments it is worth noticing that circular geodesies 
of test particles become null when cr = 1/4, for both, the static and stationary 
cylinder [10], but this corresponds to m = 1/4 in the first case and m = 1/8 
in the van Stockum case. 

On the other hand, we know of at least one example where the angular 
momentum of the source diminishes the effective active gravitational mass of 
the source: the Kerr metric. Indeed, in this later case the equation for radial 
geodesic on the symmetry axis reads: 

r = -M— — (42) 

where the dot stands for differentiation with respect to the proper time, r is 
the radial Boycr-Lindquist coordinate and a and M are the two parameters 
of the Kerr metric. It is evident from (42) that the angular momentum of 
the source a diminishes the attraction force on the test particle, as compared 
with the static case. 

All this having been said, we do not know at this point whether or not 
some kind of mechanism similar to (or different from) the one existing in the 
Kerr spacetime is acting in the van Stockum case, and more important we 
do not know what the physics underlying such mechanism, might be. 
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4 An interior solution 



In order to illustrate the use of all the equations above, let us present a simple 
model which we believe is new. As we shall see, even in a very restrictive 
situation, as the one presented below, a purely analytical treatment of the 
problem is not suitable, and one should appeal to numerical methods, for a 
complete description of the solution. Nevertheless, some information about 
the properties of the model may be obtained from the analytical approach. 
Thus, let us assume for simplicity the equations of state 

P, = afi, P, = P^ = 0, (43) 

with a = constant. When a = (43) reduces to the van Stockum solution 
[9] for a cylindrically symmetric rotating dust which is given by 

/ = k^ur^, 7 = -a;V, Z = r^(l-a;V), (44) 

«;// = 4u;V''', (45) 

where a; is a constant. 

Substituting (43) into (12) we have for 7^ 

«7' = y, (46) 

and integrating we obtain 

e"^ = C2/, (47) 

where C2 is an integration constant. Considering the regularity conditions 
(19) we can rewrite (47) as 

= /. (48) 
The field equations (6-10) with (43) reduce to 

2DD'^' + f'l' + k'^ = 0, (49) 
4(1 + a)fDD" + a \lfD^-{' - 2D^f" + 2DD' f - 3f{f'l' + k'^)] = 0. (50) 

Substituting (14), (37) and (46) into (49) and considering (11) we have 

2(1 + a)ffD' - aD{f' - Ac') = 0; (51) 
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while substituting (46) and (49) into (50) we obtain 

2(1 + a)fD" + (1 - a)Dff" - Df^ + (3 + a)ffD' = 0. (52) 
Now from (51) and (52) we have 



(/'' + 4o;c') [2(1 + a)ff" - (2 + a)f'^ + Aac' 



= 0. 



(53) 



If a = then (53) reduces to /' = 0, giving the van Stockum's solution (44). 
However if a 7^ and since 



+ 4ac' ^ 0, 



considering the transformation 



(53) reduces to 



/_ il + a)x' 
f {2 + a)x - Aac^' 



Integrating (56) we obtain 



/ = C3 (2 + a)x - 1a(? 



(l+a)/(2+a) 



(54) 

(55) 
(56) 

(57) 



where C3 is an integration constant. Prom the regularity conditions (19) and 
(22) we can rewrite (57), 



2 + q; \ (l+«)/(2+a) 



(58) 



or with the aid of (46) it becomes 



(2 + a)a 2 /2 
-/ 7 



-1 (l+a)/(2+a) 



4c2 



(59) 



Prom (59) we see that when a = it reduces to the van Stockum solution 
(44). 

Peeding back (55) into (58) we could in principle integrate to obtain 
/ = /(r), unfortunately though, the integral can only be expressed in terms 
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of hypergeometric functions, therefore further specification of the solution is 
only possible numerically. 

However, some general properties of the solution may be found without 
performing the integration of (58). 

Thus, the matching of (46) or (51) using (32-36) produces 

showing that a can have positive as well as negative values for different 
parameters that describe the source. 

Substituting (46), (49-51) and (58) into (6) we obtain 



2c 



(1 + a)/ 



(61) 



At the axis r = (61) becomes 



^,^{^)\ (62) 

If we fix the stationary parameter c we observe from (62) that if a > 0, 
or > 0, the density ji diminishes as compared to the van Stockum dust 
solution; however if o; < 0, or < 0, then /i increases as compared to the 
same solution. 

The model cannot be specified further until (58) is integrated and the 
full set of boundary conditions are used. However, as mentioned before this 
requires the introduction of numerical procedures, which is out of the scope 
of the paper. 



5 The electric and the magnetic parts of Weyl 
tensor 

The electric and magnetic parts of Weyl tensor, Ea/B and Hap, respectively, 
are formed from the Weyl tensor Cap-yS and its dual by contraction 

with the four velocity vector V^, thus 

Ea/3 = Ca-y/35VW (63) 
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where rjap-ys = +1 or — 1 for a, /3, 7, 5 in even or odd order, respectively, and 
otherwise, and with 

= ^, (65) 
= De'^. (66) 

Then it follows that there are only three non-vanishing electric components, 
namely 

E^, = (67) 

E,, = (68) 

^33 = (69) 

however they are not independent, since, by virtue of (85) they satisfy the 
relation 

Ell + E22 — ~^^Es3. (70) 

Thus, there are two independent electric components of the Weyl tensor. 

On the other hand it can be easily checked that there is only one inde- 
pendent magnetic component, namely 

H12 — — ^(^0223 — -^C'0202^ • (71) 

Accordingly, if we demand our spacetime to be purely electric {Hap = 0), 
then using (14), (71), (79) and (82) we obtain 

I'ikf - fk') = 0, (72) 

implying that either 7' = or kf — fk' = 0. In the first case it follows from 
the regularity conditions that 7 = 0. Then matching condition (35) implies 
n = 1, however this last condition implies because of (39) that cr = 0, which 
produces a Minkowski spacetime outside the cylinder. 



12 



On the other hand kf — fk' = imphes, because of regularity conditions, 
that k = 0. 

Thus purely electric cylindrically symmetric spacetimes are necessarily 
static. This is in agreement with the result by Glass [11] indicating that a 
necessary and sufficient condition for a shear-free perfect fluid to be irrota- 
tional is that the Wcyl tensor be purely electric. However, our result does 
not require the pressure to be isotropic. 

Considering the field equations (6), (7), (9) and (10) we can rewrite (67), 
(68) and (71) as 



12 



22 



1 

12 




and with (14) we have 



or using (37) 



H 



12 



ci Y 

"4 /' 



Hi2 — —- — 



c 7 
2/ 



(73) 
(74) 
(75) 

(76) 
(77) 



Since the parameter c can be shown to be proportional to the vorticity of 
the source [7], this last equation reinforces further the belief that rotation 
is one of the sources of the magnetic part of Weyl tensor [12]. However, it 
should be noticed that spacetimes are known which have vorticity but zero 
[13]. The most conspicuous example of this kind of spacetime being the 
Godel solution [14]. 



6 Conclusions 

We have deployed the equations describing the stationary (shear-free) cyhn- 
der, as well as the regularity and matching conditions and the active gravi- 
tational mass. Then a simple model has been presented to illustrate the way 
in which models can be found. 
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Next, we have calculated the Wcyl tensor components, as well as its 
electric and magnetic parts. It has been established that purely electric 
solutions are necessarily static. This result, together with a previously known 
result [1] about the nonexistence of (static) conformally flat solution which 
satisfies regularity conditions and matches smoothly to Levi-Civita spacetime 
on the boundary surface, implies the nonexistence of interior (stationary) 
conformally flat solution which satisfies regularity conditions and matches 
smoothly to Lewis spacetime on the boundary surface. 

Finally, a link between the magnetic part of the Weyl tensor and a pa- 
rameter related to the vorticity of the source has been found. 
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Appendix 



The spacetime (3) has the following non-null components of the Weyl tensor 



0101 



0202 



C'oaos — —D^e ^'''Ci2i2 - 
C'oiis 



1 

12 



f D' \ D" 



-3r + 3^/'-2-^(/'/' + 



1 

12 



/ D' \ D" 

/7"+3 /--/' y-4/- 



D' 



f 



+3r-3^/' + 4-^(/Y + fc'2) 



6 



7 



D" 

D 



1 

12 



f D' \ D" 
k-f" - ^ U— - k' U' + 2k— 



(78) 



(79) 
(80) 



14 



-3k" + 3^k'-2^{fl' + k") 



( D' \ D" 
A;7" + 3 U— -A;'J7'-4A;— 



+3^"-3^^'+4A(/'/'+0 



Cl313 — T- 



1 

12 



-3i" + 3^r-2-ij(/r + A«) 



C2323 — ~ 



1 

12 



( D' \ D" 



+31"-3^r + 4ij(/'i' + r) 



and they satisfy the relations, 

C'oiOl + Cq202 

Cq\13 + C0223 



ke^ 



0303) 



0303, 



^"1313 + C2323 — ^^(^0303, 
^C'oiOl + fC2323 — k{Co223 ~ C'oiis), 



(81) 



(82) 



(83) 



(84) 

(85) 

(86) 

(87) 
(88) 



hence we have three independent Weyl tensor components, which is one more 
than the number in the static case [1] . 
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